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ABSTRACT 

Novel  propellant  formulations  are  critical  for 
developing  Future  Force  gun  and  missile  systems  that  are 
safer,  smaller,  and  more  lethal  than  those  used  presently. 
We  employ  our  US  Army  Research  Laboratory 
combustion  model,  CYCLOPS,  to  predict  the  burning 
rates  of  nitrocellulose  with  and  without  HN3,  a  stable 
pyrolysis  product  of  high-nitrogen  compound  5- 
aminotetrazole  (5-AT).  Our  model  employs  a  detailed 
chemical  kinetic  mechanism  containing  368  chemical 
reactions  and  59  species.  We  test  its  HN3  chemistry  by 
studying  both  neat  and  HN3-doped  flames  experimentally 
and  with  the  PREMIX  flame  code.  We  measure  the  major 
and  radical  species  concentrations  by  either  molecular 
beam  mass  spectrometry,  laser  spectroscopy,  or  both,  and 
compare  them  to  those  we  model  with  PREMIX  using  the 
detailed  chemical  mechanism.  The  model  predicts  well 
the  concentration  of  the  postflame  species,  including  the 
OH  concentration,  which  decreases  by  about  12  %  with 
the  addition  of  1.1%  HN3.  However,  the  model  does  not 
predict  adequately  the  shapes  of  the  HN3  and  NH  profiles 
near  the  burner  surface.  Our  rate  and  sensitivity  analyses 
reveal  that  the  rate  expressions  of  reactions 
HN3+0H=N3+H20  and  HN3+NH=NH2+N3  should  be 
lower  by  a  factor  of  3  and  4,  respectively,  than  what  is 
reported  in  the  literature.  Our  CYCLOPS  calculations 
with  mechanism  updates  show  that  the  addition  of  HN3  to 
nitrocellulose  enhances  its  burning  rate  significantly  over 
the  10  to  300-MPa  pressure  range;  a  factor  of 
approximately  five  at  10  MPa. 

1.  INTRODUCTION 

The  ability  to  either  increase  or  decrease  a 
propellant’s  burning  rate  by  adding  selective  modifiers, 
without  otherwise  impairing  its  properties,  is  of  great 
interest  and  a  prerequisite  for  designing  novel  propellant 
formulations  with  tailored  energy  release.  One  use  of 
interest  is  to  achieve  substantial  burning  rate  differentials 
between  otherwise  similar,  compatible  propellant 
formulations  for  co-layered  propellants  in  gun 
formulations  (Leveritt  et.  ah,  2006).  Another  is  to  achieve 
very  high  burning  rates  for  rocket  applications  (Neidert, 
J.,  private  communication).  The  high-nitrogen  compound 
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triamino-guanidinium  azotetrazolate  (TAGZT)  is  one 
example  of  a  propellant  modifier.  It  has  the  strongest 
burning  rate  augmentation  effects  known  for  HMX,  RDX, 
and  nitrate  ester-based  propellants  (Flanagan,  1984; 
Walsh  and  Knott,  2004;  Conner  and  Anderson,  2008). 
The  action  of  this  burning  rate  modifier  may  be  associated 
with  small,  gas-phase  intermediate  species  such  as  NH3 
and  N2H4  that  are  produced  in  its  decomposition.  These 
species  may  be  kinetically  active  by  affecting  the  radical 
pool  growth  rate  of  the  gas-phase  flame  thereby  affecting 
the  burning  rate. 

In  this  paper,  we  report  the  results  of  our  two- 
phase,  combustion  model  CYCLOPS  for  mixtures  of 
nitrocellulose  with  HN3.  The  model  employs  a  detailed 
chemical  kinetic  mechanism  that  consists  of  over  350 
reactions  and  nearly  60  species.  We  test  key  components 
of  this  mechanism  for  accuracy  in  the  reaction’s  rate 
constants  over  a  wide  temperature  range  by  performing 
both  experimental  and  flame  code  modeling  of  neat  and 
HN3-doped  flames.  We  measure  the  species 
concentrations  by  either  molecular  beam-mass 
spectrometry  (MB/MS),  laser-induced  fluorescence  (LIF), 
or  both,  throughout  the  flames,  and  compare  them  to 
those  we  predict  with  the  PREMIX  flame  code  employing 
the  detailed  chemical  mechanism.  Rate  and  sensitivity 
analysis  reveal  key  species  and  chemical  reactions 
important  in  the  mechanism.  We  present  these  reactions 
and  the  burn-rate  effects  of  HN3  on  nitrocellulose. 

2.  EXPERIMENTAL  APPARATUS  AND 
MODELING 

Our  previous  publications  contain  the  details  of 
our  burner  apparatus,  flame  modeling,  and  CYCLOPS 
modeling  (Grams  and  Sausa,  2008;  Venizelos  and  Sausa, 
2000;  and  Miller  and  Anderson,  2004).  Briefly,  our 
burner  system  is  equipped  for  MB-MS,  LIF,  and 
thermocouple  temperature  measurements.  A  60-mm  flat 
burner  (McKenna)  mounted  within  a  cylindrical,  vacuum 
chamber  supports  the  neat  and  HN3-doped  flames  at  4 
kPa.  Vertical  translation  of  the  burner  allows  full  optical 
access  to  the  flame  through  the  chamber  view  ports.  We 
use  gas  flow  rates  of  0.86,  0.37,  and  0.71  slm  for  H2,  02, 
and  Ar,  respectively,  for  the  neat  flame,  and  0.86,  0.37, 
0.69,  and  0.02  for  H2,  02,  Ar,  and  HN3,  respectively,  for 
the  doped  flame.  The  HN3  flow  rate  corresponds  to  about 
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1.1  %  of  the  total  flow  rate.  MKS  mass  flow  controllers 
meter  the  gases,  whose  flow  rates  are  cross-checked  with 
a  GCA  wet  test  meter. 

We  model  the  flames  with  the  PREMIX  flame 
code  (ver.  2.55)  using  the  CHEMKIN  mechanism 
interpreter  (ver.  3.6)  and  chemical  kinetics  library  (ver. 
4.9)  (Kee  et  al.,  1985,  1991,  1994).  Our  calculations 
include  both  thermal  diffusion  and  multitransport  package 
options,  and  use  either  the  measured  temperature  profile 
as  input  or  solve  the  energy  equation.  The  relevant 
reactions  for  HN3  consumption  and  intermediate  and  final 
species’  formation  are  from  a  literature  review.  We  use 
the  Sandia  National  Laboratories  (SNL)  Thermodynamic 
and  Transport  Databases  for  our  calculations  (Kee  et  al., 
1985,  1985).  The  N3  and  HN3  values,  which  are  not 
included  in  the  SNL  Databases,  are  those  of  Burcat  and 
Ruscic  (Burcat  and  Ruscic,  2005). 

The  CYCLOPS  model  employs  a  detailed 
chemical  mechanism  containing  59  species  and  368 
chemical  reactions.  The  bulk  of  this  mechanism  models 
successfully  the  burning  rates  of  a  propellant  consisting 
mainly  of  nitrocellulose  (Miller  and  Anderson,  2005).  The 
present  mechanism  includes  changes  that  pertain  mainly 
to  HN3  decomposition  reactions,  and  NHX  reactions  or 
species.  These  modifications  only  modestly  change  our 
predicted  burning  rates  for  neat  nitrocellulose. 

Our  propellant  modeling  relies  on  a  pyrolysis 
law  equation  by  Zenin  and  coworkers  (Zenin  et  al.,  1995). 
Its  equation  is  as  follows:  m  (g)  =m0e  Es/RT,  where  m  (g)  is 
the  mass-flux,  T  is  the  surface  temperature,  R  is  the  gas 
constant,  m0=  1.8xl03  g/cm2  s,  and  Es  =  9935  cal/mole 
(Zenin,  1995;  Miller  and  Anderson,  2005).  We  use  the 
above  equation  for  a  given  mass  flux  and  a  set  of  surface 
products  (nascent  gas  phase  products)  to  compute  the  heat 
feedback  requirement  to  gasify  the  propellant;  the  code 
solves  the  corresponding  gas-phase,  flame  problem  and 
computes  the  heat  feedback  rate  from  the  resulting  flame 
structure  (in  particular,  the  near-surface  properties, 
especially  the  temperature  gradient  in  the  gas  phase  at  the 
surface).  The  code  performs  iterations  on  the  mass  flux, 
driven  by  an  algorithm  to  cause  the  condensed  heat  flux 
requirement  to  converge  to  the  gas  heat  flux  supplied. 
When  the  problem  converges,  we  easily  relate  the  mass 
flux  to  the  burning  rate  by  means  of  the  propellant 
density. 

3.  RESULTS  AND  DISCUSSION 

Figure  1  shows  the  flame  temperatures  of  both  neat 
and  HN3-doped  H2/02/Ar  flames  measured  with  a  thin- 
wire  thermocouple,  OH  LIF,  or  both.  The  temperature  at 


Fig.  1.  Temperature  profiles  of  an  HVO/HN }  flame 

the  burner  surface  is  about  375  K  for  both  flames.  It 
increases  from  the  burner  surface  to  about  1 1  mm  where  it 
peaks  at  -1300  K  and  then  decreases  gradually  in  the 
postflame  region.  When  we  add  about  1%  HN3  to  the  neat 
flame,  its  temperature  drops  by  about  50K  at  15.0  mm. 
Figure  1  also  shows  the  PREMIX  temperature  profile 
calculated  by  solving  the  energy  equation.  Overall,  it  is 
similar  to  the  experimental  profiles  except  in  the 
postflame  region.  The  calculated  and  experimental 
profiles  diverge  above  15  mm;  the  calculated  one 
increases,  whereas  the  experimental  one  gradually 
decreases.  This  divergence  suggests  that  the  model  does 
not  properly  account  for  the  gas-phase  heat  losses  to  the 
surroundings  when  solving  the  energy  equation.  Thus,  our 
PREMIX  calculations  use  the  experimental  temperature 
profiles  as  input  for  the  species  concentrations. 

Figure  2  shows  the  experimental  and  calculated  H2, 
02,  H20  species  profiles  for  the  doped  flame  (-1%  HN3). 
The  concentrations  of  H2,  02,  and  H20  are  similar  to 
those  in  the  neat  flame  as  expected.  Overall,  the  model 
profiles  predict  very  well  the  experimental  profiles 
throughout  the  flame.  02  reacts  completely  with  H2  to 
form  H20;  however,  H2  does  not  react  completely. 
Approximately  one  third  of  the  initial  H2  is  still  present  at 
15  mm. 

Figure  3  shows  the  experimental  and  modeled  HN3 
profiles  near  the  burner  surface  (3a,  top),  along  with  the 
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Fig.  2.  Major  species  profiles  of  an  H2/OyHN3  flame 

results  of  our  rate  analysis  (3b,  bottom).  As  can  bee  seen 
in  figure  3a,  the  model  predicts  HN3  to  react  more  rapidly 
than  it  does.  Sensitivity  analyses  of  the  HN3,  NH,  and  NO 
species  reveal  that  the  rates  of  the  following  reactions 
HN3+0H=N3+H20  (Rl)  and  HN3+NH=N3+NH2  (R2) 
influence  their  concentrations,  and  that  their  rate 
expression  should  be  lower  than  what  is  reported  by 
LeBras  and  Combourieu  and  Li  and  coworkers  (LeBras 
and  Combourieu,  1973  and  Li  et  al.,  2006).  When  we 
lower  these  reaction’s  rate  expressions  by  a  factor  of  4 
and  3,  respectively,  in  a  modified  mechanism,  the  model 
better  predicts  the  experimental  data,  including  our  NH 
data  (see  Grams  and  Sausa,  2008).  Our  experimental  data 
suggests  that  we  could  lower  the  rate  expression  of  these 
reactions  further;  however,  this  would  reduce  the 
predicted  NO  concentration  below  our  experimental 
uncertainty,  as  will  be  shown  in  the  latter  portion  of  this 
paper. 

Figure  3b  reveals  that  reactions  HN3+0H=N3+H20, 
HN3+NH=N3+NH2,  and  HN3+H=N2+NH2  account  for  the 
disappearance  of  HN3  close  to  the  burner  surface,  where 
about  90%  of  HN3  reacts  with  either  OH  or  NH.  The  heat 
release  of  reaction  HN3+0H=N3+H20  is  about  30%  that 
of  H+02+M=H02+M,  the  reaction  with  the  largest  heat 
release  in  this  system  (2.86  cal/cm3-sec). 

Figure  4  shows  both  the  N2  experimental  and 
calculated  profiles  for  the  HN3-doped  flame.  Both  models 
predict  the  N2  concentration  very  well  throughout  the 
flame.  The  profiles  show  that  the  N2  concentration 
gradually  increases  from  the  burner  surface  to  about  7 
mm,  and  then  levels  off  in  the  postflame  region.  The 
model  with  the  modified  mechanism  predicts  an  XN2/XAr 


Fig.  3.  HN3  modeled  and  experimental  profiles  (top,  3a) 
with  NH  rate  of  production  and  consumption  (bottom,  3b). 


Fig.  4.  N2  modeled  and  experimental  profiles 
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value  of  0.0442  at  19.6  mm.  This  value  is  in  excellent 
agreement  with  the  experimental  value  of  0.044.  A 
NASA-Glenn  equilibrium  calculation  yields  an  XN2/XAr 
value  of  0.0452  at  the  adiabatic  temperature,  T=2543  K. 
These  results  show  that  the  postflame  N2  concentration 
reaches  its  adiabatic  equilibrium  concentration  at 
approximately  1300  K,  and  its  concentration  is  mostly 
independent  of  temperature  in  the  temperature  range  1300 
to  2500  K. 

Presented  in  figure  5  are  both  the  NO  modeled  and 
experimental  profiles  for  the  doped  flame.  The  XN0/XAr 
experimental  value  increases  by  about  factor  of  two  from 
burner  surface  to  7  mm,  and  then  levels  off  mostly 
thereafter.  The  model  with  the  modified  mechanism 
predicts  a  value  of  0.0038  near  13  mm.  This  value  agrees 
well  with  the  experimental  value  of  0.0042  at  the  same 
location,  but  not  as  well  as  the  value  from  the  original 
mechanism.  We  employ  the  modified  mechanism,  instead 
of  the  original  mechanism  because  it  predicts  better  the 
HN3  experimental  profile,  as  shown  in  figure  3a,  and  the 
NH  profile,  not  shown  (Grams  and  Sausa,  2008).  Also,  it 
predicts  the  NO  postflame  concentration  fairly  well, 
within  the  limits  of  our  experimental  uncertainty.  A  5% 
increase  in  the  inputted  temperature  profile,  the 
uncertainty  in  our  temperature  measurement,  decreases 
the  XNO/XAr  value  by  about  2%,  whereas  a  5%  decrease  in 
the  temperature  profile  increases  the  Xno/X^  value  by 
about  2%.  These  results  reveal  that  the  PREMIX 
postflame  NO  concentrations  depend  mostly  on  the 
chemical  mechanism. 

Figure  6  shows  the  OH  experimental  and 
modeled  profiles  for  both  the  neat  and  HN3-doped  flames. 
The  experimental  profile  of  the  neat  flame  is  normalized 
to  the  modeled  profile  and  serves  as  a  basis  for  comparing 
the  modeled  OH  concentrations  to  those  determined 
experimentally  in  the  doped  flame.  The  modeled  profiles 
are  similar  to  the  experimental  profiles;  the  OH 
concentration  decreases  from  the  burner  surface  to  about 
3  mm,  increases  from  3  to  7-12  mm,  and  then  decreases  in 
the  postflame  region.  The  models  predict  well  the  overall 
shape  of  the  experimental  profiles,  including  their  decay. 
However,  the  modeled  profiles  are  shifted  about  2  mm 
away  from  the  burner  compared  to  the  experimental 
profiles.  At  present,  this  discrepancy  is  not  well 
understood.  Figure  6  reveals  that  the  addition  of  ~1%  HN3 
decreases  the  peak  OH  concentration  by  approximately 
12%.  This  decrease  is  predicted  very  well  by  the  model, 
which  shows  an  11.35%  decrease  in  OH  concentration. 

Figure  7  shows  the  experimental  burning  rates  of 
nitrocellulose  and  our  predicted  burning  rates  of 
nitrocellulose,  and  HN3  and  nitrocellulose,  employing  our 
modified  mechanism  in  our  calculations.  Our  model 
predicts  very  well  the  burning  rates  of  nitrocellulose  from 
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Fig.  5.  NO  Modeled  and  experimental  profiles 


Fig.  6.  Experimental  and  molded  OH  profiles  in  neat 
and  HN3  doped  H2/02  flame 
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about  0.5  to  350  MPa,  and  reveals  that  HN3  strongly 
increases  nitrocellulose’s  burning  rate  over  the  1  to  200- 
MPa  pressure  range,  as  much  as  seven  at  10  MPa.  A 
detailed  chemical  analysis  reveals  that  the  following 
reactions  play  a  prominent  role  in  this  augmentation 
process: 

R3:  HN3  ->NH  +  N2 
R4:  NH  +  NO  — >  H  +  N20 
R5:  NN3  +  H  ->  N3  +  H2 
R6:  HN3  +  OH  ->  N3  +  H20 
R7:  N3  +  NO  — >  N20  +  N2 
R8:  N20  +  M— >-N2  +  0  +  M 

HN3  produces  two  effects  leading  to  the 
increased  burning  rate:  one,  a  more  rapid  conversion  of 
NO  to  final  product  N2,  and  two,  some  contribution  from 
chain  branching  via  the  O-atom  radical  produced  in  R8 
which  speeds  the  kinetics.  The  relative  proximity  of  the 
heat  released  by  the  NO  conversion  to  the  surface  may  be 
the  primary  reason  for  the  very  strong,  burning  rate 
augmentation.  In  fact,  the  NO  conversion  is  predicted  to 
become  so  fast  that  the  dark  zone  that  occurs  for  pure 
nitrocellulose  at  10  atm  due  to  NO  formation  collapses 
when  HN3  is  added.  That  is,  the  first  and  second  stage 
flames  coalesce.  The  relative  proximity  of  the  heat 
released  by  the  NO  conversion  to  the  surface  as  compared 
to  NC,  where  the  dark  zone  causes  a  large  separation,  may 
be  the  primary  reason  for  the  very  strong  burning  rate 
augmentation  effect. 

CONCLUSION 

The  design  of  novel  propellant  formulations  with 
tailored  energy  release  is  important  for  Future  Force 
weapons  systems.  We  have  calculated  the  burning  rates  of 
nitrocellulose  with  and  without  F1N3  using  our  CYCLOPS 
combustion  model  with  a  detailed  chemical  kinetic 
mechanism  composed  of  368  chemical  reactions  and  59 
species.  We  tested  this  mechanism’s  FIN3  chemistry  by 
performing  a  combined  experimental  and  modeling  study 
of  both  neat  and  FlN3-doped  Fl2/02  flames.  Our  results 
show  that  the  rates  of  reactions  FIN3+0F[=N3+F[20  and 
FIN3+NFI=NF[2+N3  should  be  lower  by  a  factor  of  3  and 
4,  respectively,  than  what  is  reported  in  the  literature.  Our 
CYCLOPS  calculations  with  updated  mechanism  show 
that  the  addition  of  FIN3  to  nitrocellulose  increases  its 
burning  rate  over  a  wide  range  of  pressures  used  for 
rockets  and  guns. 
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